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The geometric distribution of high-energy hadrons >100 GeV in shower cores measured with the 
KASCADE calorimeter is analyzed. The data are checked for sensitivity to hadronic interaction 
features and indications of new physics as discussed in the literature. The angular correlation of 
the most energetic hadrons and in particular the fraction of events with hadrons being aligned are 
quantified by means of the commonly used parameter A 4 . The analysis shows that the observed 
A 4 distribution is compatible with that predicted by simulations and is not linked to an angular 
correlation from hadronic jet production at high energy. Another parameter, describing 

distances between hadrons measured in the detector, is found to be sensitive both to the transverse 
momenta in secondary hadron production and the primary particle type. Transverse momenta in 
high-energy hadron interactions differing by a factor two or more from what is assumed in the 
standard simulations are disfavoured by the measured d^^ distribution. 

PACS numbers: 96.40.Pq,96.40.-z,13.85.-t,13.85.Tp 


I. INTRODUCTION 

High-energy hadrons in extensive air showers offer a 
unique possibility to study interaction features well be¬ 
yond the kinematic and energy region of earthbound ac¬ 
celerators. In particular, structures in hadronic shower 
cores might reflect properties of the particle production 
at an initial stage of the shower development. For ex¬ 
ample, the EAS-TOP Gollaboration investigated multi¬ 
core events recorded in the calorimeter and studied the 
cross-section of large pt jet production in proton-air col¬ 
lisions P|. In other studies, data taken by the KAS- 
GADE scintillator array and hadron calorimeter Hi 
were used to investigate different aspects of hadronic in¬ 
teraction models. The correlation of the hadronic to the 
muonic and electromagnetic shower components Hi, 
as well as features of the particle production in the very 
forward region of nucleon-air collisions were measured 
and compared to model predictions. 

Geometric structures in hadronic shower cores at ob- 
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servation level are particularly interesting, as one expects 
QGD jet production to lead to secondary hadrons being 
naturally aligned to form line shape patterns 0- Simi¬ 
lar alignment structures might result from exotic hadron 
production processes 0. Therefore, it is not surprising 
that the observation of aligned structures in a number 
of emulsion chamber experiments has initiated consider¬ 
able experimental and theoretical efforts over the last two 
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Aligned event structures were, for instance, reported 
by the PAMIR experiment 0, |I3 ■ An excess of events 
with substructure alignment above background fluctua¬ 
tions was found in the data for primary energies above 
an energy threshold of 8—10 PeV (see e.g. Bin in 
Bin B | and references therein). The most energetic 
event of ~10 PeV that was measured in an emulsion ex¬ 
periment during flights with the Goncorde airplane also 
shows alignment |l (il| . Furthermore, alignment was found 
in another individual high-energy event, estimated to 
have an energy >10 PeV, that was detected by a bal¬ 
loon borne emulsion chamber B- 

The experimental evidence of the existence of an align¬ 
ment phenomenon, i.e. events with substructures being 
located in a line occurring more frequently than expected, 
is debated controversially, including the possibility of the 
excess being due to statistical fluctuations (see e.g. B)- 

At lower energy, no excess of aligned structures in air 
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showers or direct experiments was observed. At ener¬ 
gies of a few PeV the fraction of elongated events mea¬ 
sured by PAMIR was found to agree with the expecta¬ 
tions from background fluctuations. The fraction of elon¬ 
gated events of similar energy recorded by the Kanbala 
air shower experiment M gave also no indication of a 
signihcant excess over the simulated background. Mea¬ 
surements of the balloon experiment RUNJOB in a pri¬ 
mary energy range of 0.01—0.1 PeV have shown that the 
fraction of aligned events agreed with the background 
expectation [ 13 . Also the less energetic events recorded 
during the Concorde flights showed no excess alignment. 

A search for such alignment phenomena at accelerators 
was performed at comparatively low collision energy of 
250 GeV in tt—A u interactions. No signihcant excess of 
elongated events was found in data of the CERN exper¬ 
iment NA22 

Various theoretical efforts, partly hampered by the 
limited experimental guidance, were performed to ex¬ 
plain alignment by physics mechanisms. It was pointed 
out in I 3 ] that dynamical features of standard QCD 
let production might give rise to observable alignment. 
In [lOlll^ it was emphasized that the observations would 
indicate aligned particle production in the fragmentation 
region of the collision rather than in the central region. 
A connection of alignment to the production of secon¬ 
daries with high transverse momentum was suggested e.g. 
in min The appearance of an anti-shadow scattering 
mode at high collision energies as source of alignment 
was discussed in |^. An onset of semi-hard double in¬ 
elastic diffraction was proposed in |^, where alignment 
results as projection from the rupture of a quark-gluon 
string. The existence of a new particle with a long mean 
free path, produced in a new type of interaction, was also 
proposed to explain the PAMIR observation [iM 1^ . 

Simulations based on phenomenological models were 
compared to the PAMIR data in [3123 ■ It was pointed 
out that alignment features produced, e.g., in the first in¬ 
teraction of the primary cosmic ray in air can be washed 
out by subsequent interactions during the shower evolu¬ 
tion. 

While the existence of alignment features is out of 
doubt, the main questions concerning their nature are: 
To what extent can such structures observed in the data 
be related to specific particle production and interac¬ 
tion mechanisms? Are new particles or interaction chan¬ 
nels, maybe above a certain collision energy threshold, 
required to explain the data? 

In this work, geometric structures in hadronic shower 
cores recorded by the KASCADE experiment are inves¬ 
tigated. Both, the alignment of hadrons in shower cores 
and the distance between hadrons are analyzed and com¬ 
pared to detailed Monte Carlo simulations. 

The article is organized as follows. In Section QH rel¬ 
evant features of the KASCADE detector and the ap¬ 
plied data selection are described. Results on aligned 
geometric structures are presented in Sec. nni including 
a detailed study of the background expectation due to 


random fluctuations. Another observable, characterizing 
the distance between hadrons, is introduced in Sec. CYI 
and used to analyze KASCADE data. Implications and 
comparisons to results from other experiments are dis¬ 
cussed in Sec. El 


II. MEASUREMENTS 

The KASCADE air shower experiment |2 is lo¬ 
cated at an altitude of 110 m a.s.l. (corresponding to 
1020 g cm“^ atmospheric depth) at Karlsruhe, Germany. 
Its 16x20 m^ hadron calorimeter is well suited to in¬ 
vestigate geometric structures in hadronic shower cores. 
Individual hadrons with energies from 50 GeV up to more 
than 20 TeV are measured with a spatial resolution of 
~15 cm. The primary energy covered by KASCADE in¬ 
cludes the range of 8—10 PeV claimed as threshold energy 
for alignment. A detailed description of the KASCADE 
detector performance and the standard reconstruction 
procedures for the array and calorimeter data is given 

in 31. 

For the analysis, data collected by the KASCADE ex¬ 
periment in the period from May 1998 to April 2001 have 
been used. Primary energy, direction, and core position 
of air showers with energies >0.5 PeV are determined 
with the scintillator array data. The energy, position and 
direction of individual hadrons with energies >50 GeV 
are reconstructed from the measurements performed with 
the hadron calorimeter. The following selection criteria 
are applied to the data: 

(i) The reconstructed primary energy of the air shower 
{AEq/Eq ~25%) is >1 PeV, and the primary zenith an¬ 
gle is <30° (Ad ~1°). 

(ii) The reconstructed shower core location {Axc ~1.5 m) 
is contained in the active calorimeter area with a mini¬ 
mum distance to the calorimeter boundaries of 3 m. 

(iii) At least four hadrons with energies >100 GeV are 
reconstructed {AEh/Eh ~ 20 %). 

The effect of the latter condition is small even at lower 
primary energies, reducing the data set by ~13%. 

In total, 4489 events are selected this way. After trans¬ 
formation into the shower plane, the quantities A 4 and 
^max below) are evaluated for each event. An ex¬ 
ample of a measured elongated event is given in Fig. ^ 


III. ALIGNMENT OF HADRONS 

Commonly, the degree of alignment is described by a 
parameter A 4 quantifying the angular correlation be¬ 
tween the four most-energetic particles (or particle clus¬ 
ters): 
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FIG. 1: Example of a measured hadronic shower core (hadron 
positions in the shower plane) with A 4 = 0.99 and = 
5.6 m. 9 hadrons with energies above 100 GeV are recon¬ 
structed. For the four most energetic hadrons (full symbols), 
the energies in GeV are given. The shower core position as 
reconstructed by the scintillator array is marked by a cross. 
The active calorimeter area exceeds the area plotted. 


where denotes the angle between the lines connect¬ 
ing particle k to i and j. Possible values range be¬ 
tween A 4 = — i (isotropic distribution) and A 4 = 1 (per¬ 
fect alignment). Events are usually termed “aligned” or 
“elongated” for A 4 > 0 . 8 . 

The measured A 4 distribution is displayed in Fig. |21 
Also shown are the results for primary proton and iron 
nuclei simulated with the CORSIKA (v6.0) code em¬ 
ploying the QGSJET 01 hadronic interaction model. 
The simulations are performed for a continuous primary 
energy spectrum with a spectral index —2.7. It has been 
checked explicitely that adopting other values for the pri¬ 
mary energy slope or, e.g., including a knee structure 
does not influence the calculated A 4 distribution, as will 
be also clear from the following discussion. The shower 
simulation is followed by a detailed detector simulation 
based on GEANT |^. The same reconstruction algo¬ 
rithms and selection criteria are applied to simulations 
and data. The simulation statistics is of comparable size 
to the data statistics for each primary. 

Fig.H shows that the measured A 4 distribution and 
the fraction of aligned events are well reproduced by the 
standard calculations, without the need of introducing 
new particles or non-standard interaction features. More¬ 
over, despite the much smaller initial energy per nucleon 
in case of iron showers, no significant difference between 
proton and iron induced events is observed. This already 
indicates a minor sensitivity of the observed parameter A 4 
to hadronic particle production in the first interactions. 


FIG. 2: A 4 distribution measured by KASGADE compared to 
simulation results for primary proton and iron showers. For 
clarity, simulation points are slightly displaced horizontally. 



FIG. 3: Fraction of aligned events with A 4 > 0.8 as measured 
by KASGADE versus the primary energy threshold. 

The measured fraction of aligned events (A 4 > 0.8) 
is plotted in Fig. 0 as a function of the primary energy 
threshold for selecting showers. Within the statistical 
uncertainties, no energy dependence of the fraction of 
elongated events is observed. In particular, no signifi¬ 
cant increase or threshold behaviour for higher primary 
energies is found. 

In the following, the sensitivity of the observed A 4 
distribution and the fraction of aligned event structures 
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to specific particle production features is investigated. 
Events with aligned hadronic structures could be the re¬ 
sult of jet production in hadronic interactions during the 
shower development 0- The jet pair and the hadrons 
from projectile fragmentation are naturally aligned in 
one plane and the expected geometric separation falls 
in the range measurable with KASCADE. For example, 
the projected separation from the initial direction of a 
10 TeV hadron produced at an altitude of 15 km with a 
transverse momentum pt = 5 GeV is 7.5 m at KASCADE 
observation level. 

To investigate the correlation between the angular dis¬ 
tribution of particles in hadronic interactions and geo¬ 
metric structures, a simulation was carried out in which 
the azimuth angles of the produced secondary particles 
are picked at random in the centre-of-mass system of 
the collision, keeping the transverse and longitudinal mo¬ 
menta unchanged. In this way, any angular correlation of 
produced particles is eliminated. The results of calcula¬ 
tions performed for 5 PeV primary protons are shown in 
Fig. Altogether with the standard simulations. No signif¬ 
icant change of the distribution as a whole or the fraction 
of aligned events is observed for the modified version (la¬ 
beled ) with respect to the original one. 

Next, the sensitivity of A 4 to the transverse momentum 
Pt of the secondaries is checked. A connection of high- 
Pt events and alignment was suggested, e.g., in [l^ll^ . 
In the simulation, the pt of each secondary produced in 
a hadronic interaction is artificially increased by a fac¬ 
tor two (labeled “2 • pt” in Fig. 0, keeping the original 
angular correlations unchanged. Such an increase is not 
realistic, but used here to study the sensitivity. The over¬ 
all shape of the A 4 distribution is affected very little. No 
significant change of the fraction of aligned events can be 
noted. 

This insensitivity to modifications of the hadronic par¬ 
ticle production suggests that the measured A 4 distri¬ 
bution is mostly due to a random distribution of the 
hadron azimuth angles. Thus, a simplified simulation is 
performed by artificially generating sets of four hadron 
positions in the shower plane. The azimuth angles of the 
hadrons are chosen randomly with respect to the shower 
core. The lateral distances of the hadrons are sampled 
according to the measured lateral distribution of the four 
most energetic hadrons. The KASCADE data are indeed 
reproduced, see Fig.Al 

The sensitivity of the A 4 distribution to the shape 
of the lateral distribution of the observed particles is 
checked with the simplified simulation by varying the 
assumed lateral shape, keeping the choice of hadron az¬ 
imuth angles at random. It should be noted that A 4 as 
an observable, being determined from the event topology, 
does not depend on the absolute length scale of the as¬ 
sumed lateral distribution, but only on its shape. As can 
be seen in Fig. [SJ even choices of functional forms that 
are unrealistic for air shower lateral distribution such as 
adopting a particle density that is constant for any dis¬ 
tance from the core, result only in modest changes of the 
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FIG. 4: A 4 distributions for proton showers of 5 PeV pri¬ 
mary energy without modification of the hadronic particle 
production (“standard version”), with randomly chosen az¬ 
imuth angles of secondary hadrons (“‘^random”)’ with 
transverse momenta of secondary hadrons increased by a fac¬ 
tor two (“2 -pt”). Simulation points obtained after the mod¬ 
ifications are slightly displaced horizontally. 
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FIG. 5: A 4 distributions: KASGADE data and results from 
a simplified simulation of event topologies, in which the az¬ 
imuth angles were chosen randomly with respect to the shower 
core, and the distances according to the measured hadron 
lateral distribution (see text). Simulation points are slightly 
displaced horizontally. 
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FIG. 6: A 4 distributions obtained by a simplified simulation FIG. 7: distribution measured by KASCADE compared 

of event topologies (see text or caption of Fig.|^ for different to simulation results for primary proton and iron showers for 


hypothetical shapes of the lateral distribution as indicated. 


primary energies of 1-3.2 PeV. 


A 4 distribution. In general, a flatter lateral distribution 
tends to decrease the fraction of aligned events generated 
in the random sample. However, the sensitivity of A 4 to 
different lateral shapes is minor. It seems hardly possi¬ 
ble for an air shower experiment to observe differences 
between A 4 distributions that are due to different lateral 
distributions of the measured particles only. 

Other parameters describing geometric structures sim¬ 
ilar to A 4 have been checked for sensitivity to primary 
mass or hadronic interaction characteristics |3f]|. For in¬ 
stance, all hadrons >100 GeV instead of the four most 
energetic ones are considered to calculate a generalized 
quantity A„. Also, weight factors for the hadrons are in¬ 
troduced to additionally take their reconstructed energy 
into account. No significant sensitivity to hadronic par¬ 
ticle production properties of these A-related quantities 
is found, either. 

IV. DISTANCE BETWEEN HADRONS 

Hadronic shower observables with sensitivity to inter¬ 
action features have for instance been discussed in Hi 
A specific, further observable that is based on the four 
most energetic hadrons in each event as is A 4 , is given 
by the maximum distance, between one of the four 

considered hadrons to the geometric centre of the other 
three |2(i|. Since this quantity is intimately connected 
to the hadron lateral distribution, a sensitivity both to 
hadronic interaction features such as the transverse mo¬ 
mentum of secondary hadrons and to the primary particle 
type can be expected. 

The distribution measured by KASCADE at en¬ 


ergies 1-3.2 PeV is compared to simulation results for 
primary proton and iron showers in Fig. [7| The com¬ 
parison is restricted to energies below the knee due to 
the increasing uncertainty of primary flux spectra and 
composition at higher energy. The distribution for 
primary protons is peaked at smaller values compared 
to the predictions for iron-initiated events, reflecting a 
steeper hadron lateral distribution in case of air showers 
initiated by proton primaries. 

The KASCADE data are mostly bracketed by the pri¬ 
mary proton and iron expectations. This seems reason¬ 
able, as for these primary energies a mixed composition 
is favoured from extrapolations of direct cosmic-ray mea¬ 
surements, from other air shower experiments and in par¬ 
ticular from the analysis of independent KASCADE ob¬ 
servables such as the electron and muon shower sizes . 
Only at largest values, a slight underestimation 

shows up in the predictions that may indicate shortcom¬ 
ings in the modelling of hadron interactions. 

In the next step, as before, the pt of secondary hadrons 
is artificially modified in the simulations. The results are 
compared to KASCADE data in Fig.|S| Assuming trans¬ 
verse momenta of secondary hadrons produced in high- 
energy interactions twice as large as in standard simu¬ 
lations (upper panel in Fig. O, the distributions both 
for proton and iron primaries are significantly shifted to 
larger values. In this scenario, the simulations are 
hardly able to provide a satisfactory description of the 
^max Pqj. ggy assumptfon On the primary composi¬ 

tion between proton and iron nuclei, and in particular for 
a mixed composition as favoured by other measurements, 
events with small values of observed in hadronic 
shower cores are not adequately reproduced by Simula- 
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FIG. 8: distribution measured by KASCADE (see 

Fig. □ compared to simulation results for primary proton 
and iron showers with modified transverse momentum of sec¬ 
ondary hadrons. The transverse momenta were artificially 
increased (upper panel, “2 ■ pt”) and reduced (lower panel, 
“0.5 • pt’) by a factor two. 


tions with doubled pt ■ The same conclusion holds when 
artificially reducing the pt of secondary hadrons by a fac¬ 
tor two (lower panel in Fig.ISJ. In this case, the simula¬ 
tions fail to describe the events observed with large values 
of Therefore, hypothetical transverse momenta in 

high-energy secondary hadron production that differ by 
a factor two or more from the standard assumptions are 
disfavoured by the KASCADE data. 

Another application of the d™™ parameter is related 
to investigating the primary cosmic-ray composition by 


means of air shower simulations. Due to the sensitiv¬ 
ity to the primary particle type, quantities such as 
can be used to check the internal consistency of hadronic 
interaction models. 


V. DISCUSSION 

The results of the KASCADE data analysis concern¬ 
ing the A 4 parameter can be summarized as follows. 
The measured A 4 distribution and the fraction of aligned 
events of high-energy hadrons are well reproduced by the 
simulations. Moreover, the data follow the expectations 
from randomly distributed hadron azimuth angles. No 
significant dependence on the primary energy was ob¬ 
served. Within the statistical uncertainty, no correlation 
between A 4 and hadronic interaction features such as jet 
production was found. 

These results are not necessarily in contradiction to 
possible alignment excesses in other shower observables, 
at different observation levels, or in a different primary 
energy range. However, the absence of an alignment ex¬ 
cess in the present analysis might be used to constrain 
attempted alignment explanations. 

The KASCADE data can be compared to the A 4 distri¬ 
bution measured by the PAMIR detector. PAMIR [^[in| 
is an emulsion chamber experiment located at an alti¬ 
tude of 4360 m a.s.l. (corresponding to 600 g cm“^ atmo¬ 
spheric depth), measuring air shower hadrons and pho¬ 
tons with a particle energy threshold of a few TeV. 

In Fig. El the A 4 distributions obtained from KAS¬ 
CADE and PAMIR data Pin are compared for two 
different primary energies. In the energy range of a few 
PeV (upper panel of the Figure), i.e. below the thresh¬ 
old energy of the alignment excess claimed from PAMIR 
data, both distributions and the fraction of aligned events 
agree well to each other. This is a quite significant 
result, since in the PAMIR experiment, compared to 
KASCADE, a completely different detection technique 
is applied, the mass overburden is considerably smaller 
(AA„ ~ 420 g cm“^), and the energy threshold of the 
registered shower particles is higher by more than an 
order of magnitude. A possible interpretation of this 
concordance is that both A 4 distributions are generated 
mostly from a random distribution. Also possible dif¬ 
ferences in the lateral distributions of the particles mea¬ 
sured by KASCADE and PAMIR seem too small to sig¬ 
nificantly affect the A 4 distributions. 

This interpretation is supported by comparing the frac¬ 
tion of aligned events measured by KASCADE F(A 4 > 
0.8) ~ (4±1)% with the values derived from RUNJOB 
of (3±1)% [i^ and from NA22 data (~ 6 %) 1^ . 

Each result was found to agree to the respective back¬ 
ground fluctuation, and despite the completely different 
measurement conditions, the fractions of aligned events 
are comparable to each other. 

For a higher primary energy range, based on A 4 dis¬ 
tributions obtained from PAMIR data such as shown in 
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FIG. 9: A4 distributions: KASCADE and PAMIR [H El 
data. In the lower panel, the highest primary energies are 
selected. The two distributions displayed for PAMIR in the 
lower panel refer to different reconstructions from the same 
PAMIR data set El El PAMIR data are slightly displaced 
horizontally. 
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Fig. m (lower panel), an alignment excess with a thresh¬ 
old energy of 8 —10 PeV was deduced 0E1 Displayed 
in the Figure are two distributions obtained with differ¬ 
ent reconstruction methods from the same PAMIR data 
set. Included in the graph are the KASCADE results for 
primary energies above 8 PeV. Despite the completely 
different observation conditions, the distributions mea¬ 
sured by both experiments agree well to each other. The 
fraction of events with A 4 > 0.8 is larger in the PAMIR 
data compared to KASCADE. However, the statistical 
significance for the fraction of aligned events to differ be¬ 
tween the two experiments is small, amounting to ~1.5 
standard deviations only. Thus, also in this primary en¬ 
ergy range the lack of significant differences between the 
A 4 distributions, including the fraction of aligned events, 
suggests a random component as a dominant source of 
the A 4 values observed by both, PAMIR and KASCADE. 


The quantity d™™, describing distances between the 
hadrons as measured by the KASCADE calorimeter, is 
sensitive to the transverse momenta in secondary hadron 
production and to the primary particle type. Transverse 
momenta in high-energy hadron interactions differing by 
a factor two or more from the standard values are dis¬ 
favoured by the measurements. 
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